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Abstract. In a Relativistic Diffusion Model (RDM), the evolution of net- 
proton rapidity spectra with ^Jsnn in heavy systems is proposed as an indica- 
tor for local equilibration and longitudinal expansion. The broad midrapidity 
valley recently discovered at RHIC in central Au + Au collisions at J~s~n~n= 200 
GeV suggests rapid local equilibration which is most likely due to deconfinc- 
ment, and fast longitudinal expansion. Rapidity spectra of produced charged 
hadrons in d + Au and Au + Au systems at RHIC energies and their central- 
ity dependence are well described in a three-sources RDM. In central collisions, 
about 19% of the produced particles are in the equilibrated midrapidity region 
for d + Au. 
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1. Introduction 

The production and identification of a transient quark-gluon plasma in thermal 
equilibrium is of basic importance in relativistic heavy-ion physics. In this contri- 
bution I propose nonequilibrium-statistical methods to investigate analytically the 
gradual thermalization occuring in the course of stopping and particle production 
at the highest available energies. The approach is tailored to identify the fraction of 
net baryons, as well as of produced charged hadrons that attain local thermal equi- 
librium from their distribution functions in rapidity or pseudorapidity. This yields 
indirect evidence for the extent and system-size dependence of a locally equilibrated 
parton plasma. 

Net-baryon rapidity distributions are sensitive indicators for local equilibration 
and deconfinement in relativistic heavy- ion collisions [EE]- At AGS, SPS and 
RHIC energies these are clearly nonequilibrium distributions even in central colli- 
sions [Ej, and also the longitudinal distributions of produced particles are not fully 
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thcrmalized in heavy symmetric [|3] as well as in light asymmetric systems [0]. 

To account for the nonequilibrium behavior of the system, the evolution of net- 
proton rapidity spectra with incident energy ^/snn= 4.9 to 200 GeV is first studied 
analytically in a Relativistic Diffusion Model (RDM, chapter 2), and compared to 
AGS SPS [HI and RHIC [HUj central collision data, chapter 3. In addition to the 
nonequilibrium-statistical evolution as described in the RDM, collective longitudinal 
expansion is considered through a comparison of the analytical RDM-solutions with 
the data. It turns out that a three-sources RDM with a small, but sizeable (14%) 
midrapidity source is required from the RHIC net-proton data. 

The midrapidity source is more important for produced charged hadrons. In 
central d + Au collisions at RHIC energies, almost 20 % of the produced particles 
are in the equilibrium source, chapter 4. Here the centrality dependence is also 
investigated in detail. The fraction of produced particles in the midrapidity source 
tends to be larger in central Au + Au collsions [|3J- The conclusions are drawn in 
chapter 5. 



2. Relativistic Diffusion Model 



The simplified model ([G], and references therein) as discussed in this contribution 
is based on a linear Fokkcr-Planck equation (FPE) for the components k of the 
distribution function R(y,t) in rapidity space, y — 0.5 • ln{(E + p)/(E — p)), 
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Veq) ■ Rk(V,t) 



d 2 y 



D k v -R k (y,t) 



(1) 



The diagonal components Dy of the diffusion tensor contain the microscopic physics 
in the respective projectile-like (k=l), target-like (k=2) and central (k=3) regions. 
They account for the broadening of the distribution functions through interactions 
and particle creations. In the present investigation the off-diagonal terms of the 
diffusion tensor are assumed to be zero. The rapidity relaxation time t v determines 
the speed of the statistical equilibration in y-space. 

As time goes to infinity, the mean values of the solutions of Eqs. approach 
the equilibrium value y eq . This value is zero for symmetric systems, but for asym- 
metric systems like d + Au it deviates from zero. We determine it from energy- and 
momentum conservation [|SJ H>] in the system of target- and projectile- like partici- 
pants and hence, it depends on impact parameter. This dependence is decisive for a 
detailed description of the measured rapidity distributions in asymmetric systems: 



y eg (b) = 1/2 • In 



< mf(b) > exp(y max )+ < m%(b) > exp(-y max ) 



< m%{b) > exp(y max )+ < mj(b) > exp{- 



(2) 



with the beam rapidities y^ = ±y max , the transverse masses < wf 2 (&) >= 

\flrn\ 2$)+ < pt > 2 ), and masses mi^ib) of the respective participants that de- 
pend on the impact parameter b. The average numbers of participants N\ t 2(b) in 
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the incident nuclei are calculated from the geometrical overlap, or from Glauber 
theory. 

The RDM describes the drift towards y eq in a statistical sense. Whether the 
mean values of the distribution functions R± and i?2 actually attain y eq depends on 
the interaction time Tint (the time the system interacts strongly, or the integration 
time of O). It can be determined from dynamical models or from paramctrizations 
of two-particle correlation measurements. For central Au + Au at 200 A GeV, this 
yields about Ti nt ~ 10/m/c, which is too short for R\ and R2 in order to reach 
equilibrium. Note, however, that this does not apply to R3 which is born near local 
equilibrium at short times, and then spreads in time through interactions with other 
particles at nearly the same rapidity. 

The analytical diffusion model is consistent with, and complementary to parton 
cascade models where stopping involves large sudden jumps in rapidity from hard 
scatterings (eg. [0), because even hard partons can participate significantly in 
equilibration processes, as is evidenced by the high-py suppression found in Au + 
Au at RHIC. 

Nonlinear effects are not considered here. Their possible role in the context of 
relativistic heavy- ion collisions has been discussed elsewhere ([Ej, and references 
therein). These account to some extent for the collective expansion of the system 
in y— space, which is not included a priori in a statistical treatment. In the linear 
model, the expansion is treated through effective diffusion coefficients Dy** that 
are larger than the theoretical values calculated from the dissipation-fluctuation 
theorem that normally relates D y and r y to each other [[I]. This relation has been 
derived from the requirement that the stationary solution of Eq. Q is equated with 
a Gaussian approximation to the thermal equilibrium distribution in y-space. At 
fixed incident energy, this weak-coupling rapidity diffusion coefficient turns out to 
be proportional to the equilibrium temperature T as in the analysis of Brownian 
motion (Einstein relation) 

D y <x -. (3) 
T y 

One then deduces the collective expansion velocities from a comparison between 
data (yielding D^') and theoretical result (D y ), as described in more detail in [ 

na. 

The FPE can be solved analytically in the linear case with constant Dy. For 
net-baryon rapidity distributions, the initial conditions arc 8- functions at the beam 
rapidities yi, = ±y max . However, it has been shown that in addition there exists a 
central (k=3, equilibrium) source at RHIC energies which accounts for about 14% 
of the net-proton yield in Au + Au collisions at 200 AGeV [ '1 , and is most likely 
related to deconfinement. For produced particles the central source turns out to be 
even more important. 
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3. Net-proton rapidity distributions 

Comparing the solutions of Eq.QJ to Au + Au central collision (5 per cent of the 
cross section) data at AGS-energies y/s^N— 4.9 GeV [IB], it turns out that due to 
the relatively long interaction time Tint and hence, the large ratio Ti nt /T y ~ 1.08, the 
system is in rapidity space very close to thermal equilibrium (dash-dotted curve), 
with longitudinal collective expansion at tr ,,=0.49 (Fig.l, upper frame) [|2J. The 
bell-shaped experimental distribution is in good agreement with the solution of 
Eq.JTJ. The distributions remain bell-shaped also at lower energies [ |S]. This 
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Fig. 1. Net-proton rapidity spectra in the two-sources Relativistic Diffusion Model 
(RDM), solid curves, compared to central Au + Au data from AGS, Pb + Pb from 
SPS (middle), and RDM-solutions for various t/r y (bottom). From [ITT]. 

situation changes at the higher SPS energy of ^snn= 17.3 GeV (Fig. 1, middle). 
Here net-proton Pb + Pb rapidity spectra corrected for hyperon feeddown [|9] show 
two pronounced peaks in central collisions, which arise from the penetration of 
the incident baryons through the system. The gradual slow-down and broadening 
is described using Eq^ as a hadronic diffusion process with subsequent collective 
expansion, ^'',,=0.75. The associated nonequilibrium solutions with expansion for 
various values of Ti n t/T v are shown in the lower frame of Fig.l. The system clearly 
does not reach the dash-dotted equilibrium solution. Hence, both nonequilibrium 
properties, and collective expansion are required to interpret the broad rapidity 
spectra seen at the SPS. 

Within the current framework, no indication for deconfinement of the incident 
baryons or other unusual processes can be deduced from the net-proton rapidity 
data at AGS and SPS energies, because the crucial midrapidity region is here too 
small. 
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This is, however, different at RH1C energies ^/sWW= 200 GeV. The RDM 
nonequilibrium solution exhibits pronounced penetration peaks with collective lon- 
gitudinal expansion v^ oll =0.93, but it fails to reproduce the BRAHMS net-proton 
data [E3| in the broad midrapidity valley (solid curve in Fig. 2, bottom): The dif- 
fusion of the incident baryons due to soft scatterings is not strong enough to explain 
the net-baryon density in the central rapidity region. The individual nonequilib- 
rium solutions i?i and R 2 are Gaussians, and if they fit the data points near y=2, 
they necessarily grossly underpredict the midrapidity yield because it is in their 
tails. This central region can onlv be reached if a fraction of the system under- 

Central Au+Au at RHIC 



62.4 GeV 




Fig. 2. Net-proton rapidity spectra in the three-sources Relativistic Diffusion 
Model (RDM), solid curves, compared to central Au + Au data from RHIC. From 

[El. 

goes a fast transition to local thermal equilibrium, dashed curve in Fig. 2, bottom. 
With collective expansion of this locally equilibrated subsystem of 22 net protons 



,J coir 



=0.93, dashed areas), the flat midrapidity BRAHMS data are well reproduced 



in an incoherent superposition of nonequilibrium and equilibrium solutions of l(T]l 
dN(y,t = T int ) 



dy 



NxRxiy, nm) + N 2 R 2 (y, T mt ) + N eq R° q c (y) 



(4) 



with the interaction time Ti nt (total integration time of the differential equation). 
The fast transition of a subsystem of N eq ~55 baryons, or 22 protons to local 
thermal equilibrium implies that in this region of rapidity space, baryon diffusion 
has been replaced by hard scatterings of the associated participant partons with 
large rapidity transfer, and subsequent thermal equilibration. 

In a schematic calculation for the lower RHIC energy of ^snn— 62.4 GeV, I 
have used the rapidity diffusion coefficient from 200 GeV to obtain the result in the 
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upper frame of Fig. 2. Two pronounced penetration peaks can be seen, together 
with a narrow midrapidity valley. The corresponding data have been taken, and are 
presently being analyzed by the BRAHMS collaboration. Once they are available, 
an adjustment of the effective diffusion coefficient and ~N eq may be required. 

4. Particle production at RHIC 

Whereas local equilibration in rapidity space occurs only for a small fraction of the 
participant baryons in central collisions at RHIC energies, the equilibrium fraction 
tends to be larger for produced hadrons. This can be inferred from recent appli- 
cations of the RDM with three sources (located at the beam rapidities, and at the 
equilibrium value) and 5-function initial conditions [00]- After conversion to 
pseudorapidity space, the calculations by Biyajima et al. for heavy systems [ |3] 
within the three-sources RDM yield good agreement with Au + Au data at both 
130 and 200 GeV per particle pair. As an example, the result for 130 GeV is shown 
in Fig. 2 in comparison with PHOBOS data; the fit at 200 A GeV is of similar 
quality. Recently pseudorapidity distributions of primary charged particles have 
become available [1.2J as functions of centrality in d + Au collisions at a nucleon- 
nucleon center-of-mass energy of 200 GeV. They are also investigated here within 
the 3-sources RDM framework. 

For produced particles, the initial conditions are not uniquely defined. Previous 
experience with the Au + Au system regarding both net baryons and produced 
hadrons [ 0] favors a three-sources approach, with 5-function initial conditions at 
the beam rapidities, supplemented by a source centered at the equilibrium value y eq . 
Physically, the particles in this source are expected to be generated mostly from 
gluon-gluon collisions since only few valence quarks arc present in the midrapidity 
region at ^/snn — 200 GeV The final width of this source corresponds to the 
local equilibrium temperature of the system which may approximately be obtained 
from analyses of particle abundance ratios, plus the broadening due to the collective 
expansion of the system. Formally, the local equilibrium distribution is a solution 
of Q with diffusion coefficient = Dy Q , and 5- function initial condition at the 
equilibrium value. 

The PHOBOS-collaboration has analyzed their minimum-bias data successfully 
using a triple-gaussian fit [EH- This is consistent with our analytical three-sources 
approach, although additional contributions to particle production have been pro- 
posed. Beyond the precise representation of the data, however, the Relativistic 
Diffusion Model offers an analytical description of the statistical equilibration dur- 
ing the collision and in particular, of the extent of the moving midrapidity source 
which is indicative of a locally equilibrated parton plasma prior to hadronization. 

With S— function initial conditions for the Au-like source (1), the d-like source 
(2) and the equilibrium source (eq), I obtain exact analytical diffusion- model so- 
lutions as an incoherent superposition of the distribution functions Rk{y,t). The 
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Fig. 3. Charged-hadron rapidity spectra in the three-sources Relativistic Dif- 
fusion Model (RDM), solid curves, compared to central Au + Au data from 
RHIC/PHOBOS. From Biyajima et al. [Ej. 



three individual distributions are Gaussians with mean values 

< J/l,2(*) >= Ueq[l - eXp(-t/T y )] =F Umax exp (~t / ' T y ) 



(5) 



for the sources (1) and (2), and y eq for the moving equilibrium source. Hence, 
all three mean values attain y e q(b) as determined from J5J for t— > oo, whereas for 
short times the mean rapidities are smaller than, but close to the Au- and d-like 
(absolute) values in the sources 1 and 2. The variances are 



<XeM = Dl^r y [\ - exp(-2t/r y 



(G) 



The charged-particle distribution in rapidity space is then obtained as incoher- 
ent superposition of nonequilibrium and local equilibrium solutions of as in Q , 
with N k replaced by N^ h . The average numbers of charged particles in the Au- and 
d-like regions are proportional to the respective numbers of participants iVi,2j 



N 



ch 



(N 1 +N 2 ) 



(7) 



with the constraint AT'jJ* 



- N^ h + N^?. Here the total number of charged 
particles in each centrality bin N*ff is determined from the data. The average 
number of charged particles in the equilibrium source Nf? is a free parameter that 
is optimized together with the variances and Ti nt /T y in a x 2 ~St of the data using 
the CERN minuit-code. 

The results for N±,2 in a geometrical overlap calculation are consistent with 
the Glauber calculations reported in [ I12| which we use in the further analysis. 
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The corresponding equilibrium values of the rapidity vary from y eq = - 0.169 for 
peripheral (80-100%) to y eq = - 0.944 for central (0-20%) collisions. They are 
negative due to the net longitudinal momentum of the participants in the laboratory 
frame, and their absolute magnitudes decrease with impact parameter since the 
number of participants decreases for more peripheral collisions. 

The result of the RDM calculation is shown in Fig. 4 (right column, middle) for 
central collisions (0-20%) of d + Au. The charged-particle yield is dominated here 
by hadrons produced from the Au-like source, but there is a sizeable equilibrium 
source (dash-dotted) that is more important than the d-like contribution. This 
thermalized source is moving since y eq has a finite negative value for d + Au, 
whereas it is at rest for symmetric systems. 




Fig. 4. Charged-hadron rapidity spectra in the three-sources Relativistic Diffusion 
Model (RDM), solid curves, compared to d + Au data from RHIC/PHOBOS for 
various centralities. From [114). 

The total yield is compared to PHOBOS data [ Es] which refer to the pseu- 
dorapidity r) = — ln[tan(0/2)] since particle identification was not available. As a 
consequence, there is a small difference to the model result [0] in y-space (y « i]) 
which is most pronounced in the midrapidity region. It is removed when the theo- 
retical result is converted to r;-space through the Jacobian 

J( V , (m)/(pr)) = coshfa) • [1 + ((m)/(p T )) 2 + sinh 2 ^)]- 1 / 2 . (8) 

Here the average mass < m > of produced charged hadrons in the central region 
is approximated by the pion mass m w , and a mean transverse momentum < px > 
= 0.4 GeV/c is used. In the Au-like region, the average mass is larger due to the 
participant protons, but since their number is small compared to the number of 
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produced charged hadrons in the d + Au system, the increase above the pion mass 
remains small, and it turns out to have a negligeable effect on the results. 

The magnitude of the equilibrium source for particle production is about 19% 
in the light and asymmetric d + Au system, whereas it tends to be larger in the 
Au + Au system 

The model calculations are compared with PHOBOS data for five centrality 
cuts [El and minimum bias [Ej in Fig- 4. The observed shift of the distributions 
towards the Au-like region in more central collisions, and the steeper slope in the 
deuteron direction as compared to the gold direction appear in the Relativistic 
Diffusion Model as a consequence of the gradual approach to equilibrium: The 
system is on its way to statistical equilibrium, but it does not fully reach it. 

Given the structure of the differential equation that is used to model the equi- 
libration, together with the initial conditions and the constraints imposed by Eqs. 
(J2J and |(7J|, there is no room for substantial modifications of this result. In partic- 
ular, changes in the impact-parameter dependence of the mean values in Eq. (JSJ 
that are not in accordance with Eq. J5J vitiate the precise agreement with the data. 

5. Conclusion 

To conclude, I have interpreted recent results for central collisions of heavy systems 
at AGS, SPS and RHIC energies in a Relativistic Diffusion Model (RDM) for multi- 
particle interactions based on the interplay of nonequilibrium and local equilibrium 
(" thermal" ) solutions. In the linear version of the model, analytical results for the 
rapidity distribution of net protons in central collisions have been obtained and 
compared to data. 

The enhancement of the diffusion in rapidity space as opposed to the expecta- 
tion from the weak-coupling dissipation-fluctuation theorem has been interpreted 
as collective expansion, and longitudinal expansion velocities for net protons have 
been determined from a comparison between RDM-results and data based on a rel- 
ativistic expression for the collective velocity, cf. f° r details. A prediction for 
net-proton rapidity distributions at ^/Sjvjv= 62.4 GeV yields a smaller midrapidity 
valley than at 200 GeV which will soon be compared with forthcoming data. 

Charged-particle production in d + Au collisions at ^/snn= 200 GeV as func- 
tion of centrality has also been investigated within the framework of my analytically 
soluble three-sources model, and compared to Au + Au. Excellent agreement with 
recent PHOBOS pseudorapidity distributions has been obtained. For central col- 
lisions, a fraction of only 19% of the produced particles arises from the locally 
equilibrated midrapidity source, whereas this figure tends to be larger in Au + Au 
collisions at the same energy. 

At y/SNN= 200 GeV, the midrapidity source for particle production must es- 
sentially have partonic origin, because the net-proton rapidity spectra for Au + Au 
show that the nucleon (and hence, valence-quark) content in this region of rapidity 
space is small Since the source is in local thermal equilibrium, it is likely that 



10 



G. Wolschin 



the generating partons form a thcrmalizcd plasma prior to hadronization. 

The d + Au results show clearly that only the midrapidity part of the distribu- 
tion function comes very close to thermal equilibrium, whereas the interaction time 
is too short for the d- and Au-like parts to attain the thermal limit. The same is 
true for the heavy Au + Au system at the same energy, but there the precise frac- 
tion of particles produced in the equilibrium source is more difficult to determine 
due to the symmetry of the problem. 
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